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The 3 + 1 formalism of Thorne and Macdonald has been used to derive the linear two-fluid equations 
for transverse waves in the plasma closed to the Schwarzschild-anti-de Sitter (SAdS) black hole. We 
reformulate the relativistic two-fluid equations to take account of gravitational effects due to the event 
horizon and negative cosmological constant and describe the set of simultaneous linear equations for the 
perturbations. Using a local approximation we investigate the one-dimensional radial propagation of 
Alfven and high frequency electromagnetic waves. We derive the dispersion relation for these waves and 
solve it for the wave number k numerically. 

Keywords: Two-fluid plasma, Alfven and high frequency electromagnetic waves, Black hole in anti-de 
Sitter space 
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O^Black holes belong to the most fascinating objects pre- 
5^ dieted by Einstein's theory of gravitation. Black holes 
bJQare still mysterious p]. In recent years the investigation 
of plasmas in the black hole environment is important be- 
cause a successful study of the waves and emissions from 
plasmas falling into a black hole will be of great value in 
aiding the observational identification of black hole candi- 
dates. The theory of general relativity and its application 
to the plasma close to the black hole horizon, however, 
have remained esoteric, and little concrete astrophysical 
impact has been felt. Within 3i? s (3 Schwarzschild radii) 
it is possible to have plasma [21 02 IH [5]. The plasma 
in the black hole environment may act as a fluid and 
Black holes greatly affect the surrounding plasma medium 
(which is highly magnetized) with their enormous grav- 
itational fields. Hence plasma physics in the vicinity of 
a black hole has become a subject of great interest in 
astrophysics. 

Sakai and Kawata have developed the linearized treat- 
ment of plasma waves using special relativistic formula- 
tion. Such an investigation of wave propagation in general 
relativistic two-fluid plasma near a black hole is impor- 
tant for an understanding of plasma processes. In order 
to investigate the behavior of plasma, a more intuitive 
mathematical description of general relativistic two-fluid 
plasma has been introduced (Buzzi, Hines, and Treumann 
(BHT)) [6] [7] based on the 3+1 membrane paradigm 
(Thorne, Price, and Macdonald (TPM)) [H EJ \M 03] near 
the horizon of the Schwarzschild black hole. In this paper 
we shall apply this method to investigate the nature of 
the two-fluid plasma (electron and positron or ion) near 
the SAdS black hole atmosphere. 

The membrane paradigm is mathematically equiva- 
lent to the standard, full general relativistic theory of 
black holes in which the black holes event horizon are 
replaced with a membrane endowed (electric charge, elec- 
trical conductivity, and finite temperature and entropy) 
so far as physics outside the event horizon is concerned. 



But the formulation of all physics in this region turns 
out to be very much simpler than it would be using the 
standard covariant approach of general relativity. In the 
TPM formulation, work connected with black holes has 
been facilitated by replacing the hole's event horizon with 
a membrane endowed with electric charge, electrical con- 
ductivity, and finite temperature and entropy. 

The study of plasma wave in the presence of strong 
gravitational fields using the 3 + 1 approach of is still 
in its early stages. Zhang [T2] [T3] has considered the 
care of ideal magneto hydrodynamics waves near a Kerr 
black hole, accreting for the effects of the holes angular 
momentum but ignoring the effects due to the black hole 
horizon. Holcomb and Tajima [14 , Holcomb [15], and 
Dettmann et al. |16j have considered some properties of 
wave propagation in a Friedmann universe. 

Actually the 3 + 1 approach was originally developed 
by Arnowitt, Deser, and Misner [17] to study the quanti- 
zation of the gravitational field. Since then, their formu- 
lation has most been applied in studying numerical rela- 
tivity [18] . TPM extended the 3 + 1 formalism to include 
electromagnetism and applied it to study electromagnetic 
effects near the Kerr black hole. 

The principal objective of this work is to solve the 
linearized two-fluid equations for the perturbations in the 
Rindler coordinate and then to determine the Alfven and 
high frequency electromagnetic waves by using the local 
approximation with appropriate boundary conditions in 
the vicinity of Schwarzschild-anti-de Sitter Black Hole. 

The solutions of black holes in Anti-de Sitter spaces 
come from the Einstein equations with a negative cos- 
mological constant. Anti-de Sitter black holes are differ- 
ent from de Sitter black holes. The difference consist- 
ing in them is due to minimum temperatures that occur 
when their sizes are of the order of the characteristic ra- 
dius of the anti-de Sitter space. For larger Anti-de Sitter 
black holes, their red-shifted temperatures measured at 
infinity are greater. This implies that such black holes 
can be in stable equilibrium with thermal radiation at a 



certain temperature. Recently the AdS/CFT correspon- 
dence suggested that String theory in anti-de Sitter space 
(AdS) equal to Conformal Field Theory (CFT) in one 
less dimension [THJ [2U1 HI] an d there is a duality between 
quantum gravity in an asymptotically de Sitter spacetime 
and a CFT living on its boundary g2 HQ H2 OH [27] . 
So our study on the Schwarzschild-anti-de Sitter black 
hole is reasonable and meaningful. 

This paper is organized as follows. In section [2] we 
summarize the 3 + 1 formulation of general relativity. 
We describe the nonlinear two-fluid equations expressing 
continuity and conservation of energy and momentum in 
section [3] The two-fluids are coupled together through 
Maxwell's equations for the electromagnetic fields. For 
zero gravitational fields these equations reduce to the cor- 
responding special relativistic expressions. In section [4] 
we restrict one-dimensional wave propagation in the ra- 
dial z direction. We linearized the two-fluid equations for 
transverse wave propagation with Rindler approximation 
in section [5j In section [5] we express the value of the un- 
perturbed quantities. We discuss the local or mean-field 
approximation used to obtain numerical solutions for the 
wave dispersion relations. We describe the dispersion re- 
lation for the transverse waves in section and give the 
numerical procedure for determining the roots of the dis- 
persion relation. The Alfven and high frequency modes 
are shown in sections [8] and [9] Finally, in section [10] we 
present our remarks. We use units G = c = kg = 1. 

2 Formulation of 3 + 1 Spacetime 

The 3 + 1 formulation of general relativity developed by 
TPM [5] E3 [TUJ ECE] is based on the concept of selecting a 
preferred set of spacelike hyper surfaces which form the 
level surfaces of a congruence of timelike curves. The 
metric for a static spherically symmetric spacetime with 
mass M and a negative cosmological constant A = — 3/£ 2 
in asymptotically anti-de Sitter spacetime has the form 

ds 2 = -f(r)dt 2 + -— dr 2 + r 2 (d0 2 + sm 2 dd(p 2 ), (1) 
f( r ) 



where the metric function f(r) is 

ti \ i 2M 7-2 



(2) 



and the coordinates are defined such that — oo < t < oo, 
r > 0, < 6 < 7T and < 4> < 2n. The function f(r) 
vanished at the zeros of the cubic equation 



r 3 + £ 2 r - 2M£ 2 = 0. 
The only real roots of this equation is 



r+ 



' : VS£ sinh[- sinh -1 (3^3-7- )]. 



(3) 



(4) 



Expanding r + in terms of M with 1/i 2 « M 2 /9, we 
obtain 

4M 2 

r + =2M(l- — + ....). (5) 



Therefore, we can write r + = 2 M rj, with 77 < 1. The 
event horizon of the SAdS black hole is smaller than the 
Schwarzschild event horizon, ru = 2M. An absolute 
three-dimensional space defined by the hyper surfaces of 
constant universal time t is described by the metric 



ds 2 



1 



dr 2 +r 2 (d9 2 + sin 2 Odip 2 ). 



(6) 



The fiducial observers (FIDO) at rest in this absolute 
space measured all local physical quantities using FIDO's 
proper time. For this FIDO use a local Cartesian co- 
ordinate system with unit basis vectors tangent to the 
coordinate lines 



3f = VW) 



d_ 

dr ' 



ld_ 

rd9' 



_L d_ 

ramd dip 



(7) 



For a spacetime viewpoint rather than a 3 + 1 split of 
spacetime, the set of orthonormal vectors also includes 
the basis vector for the time coordinate given by 



d _ 1 d 
dr a dt ' 



(8) 



where a is the lapse function (or redshift factor) defined 

by 

dr , 2M r 2 , , /9 
<r) = — = {l + ^) 1/2 - (9) 



The gravitational acceleration felt by a FIDO is given 

by [sucroiin] 

1 ,M 2r\ , , 

a=-Vlna = (_ + _) ef) 10 

while the rate of change of any scalar physical quantity or 
any three-dimensional vector or tensor, as measured by a 
FIDO, is defined by the derivative 



d fid 



(11) 



v being the velocity of a fluid as measured locally by a 
FIDO. 



3 Two-fluid Equations 

In this section we consider the plasma mixture of two 
perfect fluids either electron-positron or electron-ion and 
describe the two-fluid equations for continuity, the con- 
servation of energy and momentum, and Maxwell's equa- 
tion in 3 + 1 formalism. Here we choose these type of 
fluid component because the dispersion relation that re- 
sult from the following investigation are valid for either 
one of this two fluids since no assumption can be made 
upon the mass, number density, pressure or temperature 
of this two fluids. In the TPM 3 + 1 formulation, the 
continuity equation for each of the fluid species is 



d 

Tr:{lsn s ) + V • (a7 s n s v s ) = 0, 
dt 



(12) 



where s is 1 for electrons and 2 for positrons (or ions). For 
a perfect relativistic fluid of species s in three-dimensions, 



the energy density e s , the momentum density S s , and 
stress-energy tensor W| fc are given by 

e s = 7 2 (e s + P s v 2 ), S s = 7 2 ( £s + P s )y s , 
W* = 1 2 s (e s +P s )viv k s +P s gi k , (13) 

where v s is the fluid velocity, n s is the number density, P s 
is the pressure, and e s is the total energy density defined 

by 

e s = m s n s + P s /(l g - 1). (14) 

The gas constant j g take the value 4/3 for T — > oo and 
5/3 for T — > 0. The ion temperature profile is closely adi- 
abatic and it approaches 10 12 K near the horizon [28]. Far 
from the (event) horizon electron (positron) temperatures 
are essentially equal to the ion temperatures, but closer 
to the horizon the electrons are progressively cooled to 
about 10 8 — 10 9 K by mechanisms like multiple Compton 
scattering and synchrotron radiation. Using the conser- 
vation of entropy the equation of state can be expressed 

where D/Dt = (l/a)d/dt + v s • V. The full equation of 
state for a relativistic fluid, as measured in the fluid's rest 
frame, is as follows [2"9l IBTfl: 



iH^\im s n s /P s ) 
Lff^ 1) (im,n./P,). 



(16) 



where the H^\x) are Hankel functions. The quantities 
of Eq. (TT3"]) for the electromagnetic field take the following 
form: 

es = ^(E 2 + B 2 ), S s = i E xB, 



1 



B 2 V fc - -r{E>E k + B 3 B k ).(17) 

47T 



The conservation of energy and momentum equations are 
written, respectively, as follows [SI HI EH HI] : 



ld_ 

a dt 

a dt 



-VS s + 2a-S s , (18) 
e s a- -V- (aW,). (19) 



When the two- fluid plasma couples to the electromagnetic 
fields, Maxwell's equations take the following form: 



V B 

V ■ E 

dB 

~dt 
<9E 

~dt 



0. 

47T(T, 



(20) 
(21) 



-V x («E), (22) 
V x (aB) - 4vraJ, (23) 
where the charge and current densities are defined by 

<J = ^ Tsgs^s, J = ^7 s <j s n s v s . (24) 



Using Eq. (|13p. the energy and momentum conservation 
Eqs. ([T8]) and fp~9|) can be rewritten for each fluid species 
s in the following form 



a dt a at 



-7 s g s n s E • v s + 2^1 (e s + P s )a • v s = 0, 



r:u, - /' J ( ^§i + v s - V ) v. - v/'. 



(25) 



- 7fl g fi n s (E + v s x B) + v s 7s g s n s E • v s + -^--P, 
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+7 2 (e s + P s )[v s (v s -a) -a] = 0. 



(26) 



Although these equations are valid in a FIDO frame, for 
a = 1 they reduce to the corresponding special relativistic 
equations as given by SK [31) which are valid in a frame 
in which both fluids are at rest. The transformation from 
the FIDO frame to the commoving (fluid) frame involves 
a boost velocity, which is simply the freefall velocity onto 
the black hole, given by 



v s = (1- a 2 )" 



(27) 



Then the relativistic Lorentz factor becomes 7 boost — (1 — 

4 Radial Wave Propagation 

The two-fluid equations describing transverse and longi- 
tudinal waves can be separated by considering one dimen- 
sional wave propagating in radial z direction and intro- 
ducing the following complex variables 

v sz (z,t) = u s (z,t), v s (z,t) = v sx (z,t) + iv sy (z,t), 
B{z,t) = B x (z,t) + \B y {z,t) 1 

E(z,t) = E x (z,t)+iE y (z,t). (28) 

Then we have 

v sx B y ~v sy B x = ^(v s B*-v*B), 

v sx E y -v sy E x = ^(v s E*-v*E). (29) 

where the * denotes the complex conjugate. The conti- 
nuity equation (fT2|) and Poisson's equation ([2T]) become 



d d 

— (j s n s ) + —(a-f s n s u s )=0, (30) 



dE z 
dz 



= 47r(<7ini7i + q 2 n 2 ^ 2 ). 



(31) 



The e x and eg components of Eqs. (|2"2"j) and (|23l) give 

18B = _ i (d__ Q , £ 
a dt \dz 
dE ( d 

i— = -a ( — ~ a) B - i4nea(-y 2 n 2 v 2 - 7i"i u i) (33) 



(32) 



Differentiating Eq. ([33| with respect to t and using Eq. 
(l32l). we obtain 



d 2 



da d d 2 ,da. 



dz 2 dz dz dt 2 ^ dz ^ J 
d 

= 47rea— (n 2 72«2 - «i7i^i)- (34) 
dt 



From the e± and components of the equation 
the transverse component of the momentum conservation 
equation is 

Dv 

Ps— 1 = q s n s ^ s (E - \v s B z + \u s B) - u s v s p s a 
Dt 

( 1 dP s \ 

—V„ I q s n s j s 'E ■ v s + —-qT- J , (35) 

where E • v s = \{Ev* s + E*v s ) + E z u s , and p s is the 
total energy density defined by 

Ps = 7 2 (e s + P s ) = -f 2 {msn s + T g P s ), (36) 

with T g = -f g /{jg - !)• 

5 Linearization and Rindler 
Approximation 

In this section we linearized the two-fluid equations for the 
transverse waves using the following definitions as BHT 

u s (z,t) = uq s (z) + Su s (z,t), 

n s (z,t) = n 0s (z) + 5n s (z,t), 

P s (z,t) = Pqs(z) + 6P s (z,t), v s (z,t) = 6v s (z,t), 

p s (z,t) = p as (z) + 5p s (z,t), E(z,t) = SE(z,t), 

B z (z,t) = B (z) + SB z (z,t), B(z,t) = SB(z,t)(37) 

Here, magnetic field has been chosen to lie along the ra- 
dial e z direction. The relativistic Lorentz factor is also 
linearized such that 

7s = 70s + ^7s, where 

7o s = (1 - u^ s ) 2 , 5js = 7os u os • Su s . (38) 

The unperturbed radial velocity near the event horizon 
for each species as measured by a FIDO along e z is as- 
sumed to be the freefall velocity so that 

u 0s (z) =v s (z) = [l-a 2 (z)]i (39) 

The conservation of entropy, Eq. (1151) becomes 



6Ps = lE^lSn. 



(40) 



and from the total energy density, Eq. (f5B|) redice to 
Sps = ^(l+ 7 ° 2s7gP °0 Sn s + 2u Q sllsPvs5u s , (41) 

HQs \ POs J 



where po s = 7o s (m s no s + r s Pos)- Using Eq. [37] the trans- 
verse part of the momentum conservation equation is lin- 
earized to 



d_ d da i iaq s josnosB \ dSv. 

dz c 
aq s j 0s n 0s 



ot-U 0s — + ^7 - «o» r , 
dz at dz 



Po. 

d d da 



dt 



Poisson's equation (j3"Tj) and Eq. (|34l) are linearized to 
obtain, respectively 



d5E z 
~dz~ 



47re(n 02 702 - «oi7oi) + -iTre(j 02 Sn 2 - 7oi<5™i) 



+47re(n 02 M 27o2' 5 '"2 - noiUoi7oi^i), 



(43) 



dSv2 dSvi 



dz 2 



4irea n 02 702- 



dt 



«oi7oi- 



dt 



(44) 



The Rindler coordinate system, in which space is locally 
Cartesian, provides a good approximation to the SAdS 
metric near the event horizon. In this coordinates the 
two-fluid equations can be solved just inside and outside 
the event horizon and cannot be used to the case for ex- 
tremal black holes because there are no Rindler coordi- 
nates locally near horizons. Near the event horizon of 
SAdS black hole, choosing the coordinate 



/«(r-r+)x2( 



(r+/£) 2 + Mt/r 2 



) = M 2 , (45) 



we may write the matric ([lj approximately as 

{in) 2 



ds 



-{Kz) z dr+ 



(r+/£ + M£/r 2 + ) 2 



dz 2 +r 2 {d0 2 +sm 2 6dcp 2 ) 
(46) 

This matric becomes a Rindler spacetime with the surface 
gravity 

l/r+ M£\ 



A l 



(47) 



In the Rindler coordinates the SAdS spacetime becomes 
ds 2 = -a 2 dt 2 + dx 2 + dy 2 + dz 2 , (48) 

where 

x = r+ (6 - £ ) , y = r+ip 



z = -. (49) 

K 



The standard lapse function in Rindler coordinates be- 
comes a — zk. 

6 Radial Dependence of 
Equilibrium Parameters 

The dependence of the equilibrium field and fluid param- 
eters on the radial coordinate for transverse waves is dis- 
cussed here. The equilibrium parameters are expressed 



in terms of limiting horizon values and then for the lapse 
function and the equilibrium fields and fluid quantities 
a local (or mean-field) approximation is used. From the 
continuity equation (|30[) . we have 

r 2 ajosnosUQs = const. = r 2 _a + j + n + u + , 

where the subscript + indicate the limiting values at the 
event horizon. The freefall velocity at the horizon be- 
comes unity so that u + = 1. Since uq s = vs, 70s = 
l/a; and hence ajo s — a+7+ = 1. Also, because vg = 
(r+/r) 1 / 2 C with < C < 1- The value for C = 1 cor- 
responds to the Schwarzschild black hole. The number 
density for each species can be written as follows: 



n 0s (z) = n +s v s {z)/C 4 , 



(50) 

The equation of state Eq. (fT5|) lead to write the unper- 
turbed pressure, 



n s 



P 0s (z) = P +s (^ 



(51) 



In terms of the freefall velocity, it can be written as 

P 0s (z) = P +s v^(z)/C 16/3 . (52) 



Since Pq s = fcsno s 2o s , then with ks = 1, the temperature 
profile is 

To^T+.vp'-Viz)/?'*. (53) 

The unperturbed magnetic field is purely in the radial 
direction. It does not experience effects of spatial curva- 
ture. From the flux conservation equation V • Bo = we 
have 

r 2 B {r) = const., 

in terms of freefall velocity one can obtains the unper- 
turbed magnetic field in the form 



B (z) = B + v 4 s (z)/C 4 



(54) 



where vg = [1 — a 2 (z)] 1 ' 2 . 

In this work our main subject is how the waves behave 
in the vicinity of black hole and how the horizon affects 
these waves in anti-de sitter spacetime. Since the plasma 
is situated relatively close to the horizon, a 2 -C 1, then a 
relatively small change in distance z will make a signifi- 
cant difference to the magnitude of a. We need to choose 
a sufficiently small range in z for which the values of a 
does not vary much. So we consider thin layers in the e% 
direction, each layer with its own ao, where ao is some 
mean value of a within a particular layer. Then a more 
complete picture can be built up by considering a large 
number of layers within a chosen range of ao values. 

The local approximation imposes the restriction that 
the wavelength must be smaller in magnitude than the 
scale of the gradient of the lapse function a, i.e., 



A < 



da 
dz 



= 2r 4 



5.896C x 10 5 cm, 



or, equivalently, 
2r+ 



1.067C" 1 x lO^cm^O < C < 1, 



for a black hole of mass ~ 1M Q . 

In the local approximation for a, a ~ ao is valid 
within a particular layer. Hence, the unperturbed fields 
and fluid quantities and their derivatives, which are func- 
tions of a, take on their corresponding "mean- field" values 
for a given ao. Then the coefficients in Eqs. (|48|) . (|4"TT) . 
and (|52| are constants within each layer with respect to 
a (and therefore z as well). So it is possible to Fourier 
transform the equations with respect to z, using plane- 
wave-type solutions for the perturbations of the form ~ 
e i(kz-ut) f or eacn aQ i averp 

Using Fourier transformation, Eqs. (|51l) and (1551) become 



5E 



i47rea u;(7io27o2^2 - »-oi7oi^i) 
aok(aok — i3n) — uj 2 — k 2 



(55) 



, , . a q s j 0s n 0s B Q \ 

to aokuo s — bj + iuq s k H dv s 

POs J 

-iao (aokuos — to — iuq s k) oh = 



f>0s 



(56) 



7 Dispersion Relation and 

Numerical Solution of Modes 

The dispersion relation for the transverse electromagnetic 
wave modes may be put using Eqs. (j5"5"j) and in the 
form 



[K± {K± ± in) -uj 2 + k 2 ] = a 2 < 

U 2 2 (u-U()2K±) 1 



(u 02 K T + a w c2 ) 



ujpxiuj - U lK ± ) 
(u iK T - u — a u) c i) 

(57) 



for either electron-positron or electron-ion plasma. 

Here the local plasma frequency for transverse wave 
ojps — \/ 47re 2 7Q S 7iQ S /pos, depends on the local number 
density of electrons and also upon the local value of the 
lapse function a, and the local cyclotron frequency uj cs = 
e 2 jo s no s Bo/ pos, depends upon both the local magnetic 
field and the lapse function, with K± = aok ± i/2r+. 
The + and — denote the left (L) and right (R) modes, 
respectively. The dispersion relation for the L mode is 
obtained by taking the complex conjugate of the disper- 
sion relation for the R mode. The two modes have the 
same dispersion relation in the special relativistic case. 
The dispersion relation we have derived here for trans- 
verse waves are complicated enough and an analytical so- 
lution is cumbersome and unprofitable, even in the sim- 
plest cases for the electron-positron plasma where both 
species are assumed to have the same equilibrium param- 
eters. We solve numerically the dispersion relation in the 
form of a matrix equation as follows: 



(A - H)X = 0. 



(58) 



The eigenvalue is chosen to be the wave number k, the 
eigenvector X is given by the relevant set of perturba- 
tions, and / is the identity matrix. We need to write 



the perturbation equations in an appropriate form. We 
introduce the following set of dimensionless variables: 

kc 



w 

oj = 

j.. Su s 

0U S = 

UOs 

-DO 



k — , fc+ — , 

CJ* id* 

Sv s Sn s 

v s = , on s = 

uos n 0s 

~ SE ~ SE Z 

E= — , 5E Z = — -. 

Bo Bq 



(59) 



For an electron-positron plasma, u) p \ — ui P 2 and uj c \ = 
uici'i hence, is defined as 

lo c Alfven modes, 

high frequency modes, 

where lo p = ^Ju p ilu P 2 and u> c = ^/w c iw C 2. However, for 
the case of an electron-ion plasma, the plasma frequency 
and the cyclotron frequency are different for each fluid, 
and so the choice of cj* is a more complicated matter. For 
simplicity, we assume that 



h(^ci + ljj2 c2) 1 ^ 2 Alfven modes, 



where 



V2 



a/2 



high frequency modes, 



(2w2 +u; 2 cs ). 



The dimensionless eigenvector for the transverse set 
of equations is 



transverse 



Six 

8i>2 

SB 
SE 



(60) 



When Eqs. (|32|) and (|33|) are linearized and Fourier 
transformed, they take the forms 



IK 

O'O 



— SE 
-HoinoiSvi). 



SE 
k- 



Q,) 



SB = 0, (61) 
SB + 47re(7o2no2i5w2 

(62) 



Using Eq. ([55]), we write Eq. dHTJ), and in the 

dimensionless form: 



kSv, = ( — 







<-^cs 


UOs 


(I) 


UOs^O* 






SB-i 


(1) 


Uo s LO* 



- ^ 1 

ao 



(I) 



kSB 



"oi- 



Svi - u 02 - 



Sv 2 



ik e 
ao 



(63) 
(64) 

-SB + iCbSE 



(65) 



These are the equations in the required form to be used 
as input to Eq. ([57]). Equations ([55]). (IMl) . and (JBSJ) can 
be written as (A — kI)X = 0. We have calculated all 
the eigenvalues k of the complex matrix A to draw the 
Alfven and high frequency transverse wave modes using 
MATLAB. 






Figure 1: Top: Real part of the complex conjugate pair 
of Alfven modes for electron-positron plasma. Center: 
Imaginary part of growth mode. Bottom: Imaginary part 
of damped mode 

8 Alfven Modes 

1. Electron-positron Plasma 

For the ultrarelativistic electron-positron plasma in the 
special relativistic formulation, only one real Alfven mode 
is found to exist [31, this is because both the left and right 
circularly polarized modes were described by the same 
dispersion relation thereby leading to the same mode and 
BHT found two Alfven modes for general relativistic for- 
mulation close to the horizon of a Schwarzschild black 
hole. In our present work we have found four modes 
for electron-positron plasma two of which are shown in 
Fig. are complex conjugate pair and are similar with the 
two modes of BHT [6]. The third mode is damped sim- 
ilar to the damped mode of Fig. other two modes are 
also damping and growth modes and have a little differ- 



ence in with smaller damping rate and the fourth mode is 
equivalent to the third mode with opposite in sense and 
shows growth. The damping and growth rates of all the 




Figure 2: Top: Real part of the complex conjugate pair of 
Alfven modes for electron-ion plasma. Center: Imaginary 
part of damped mode. Bottom: Imaginary part of growth 
mode 

modes are smaller by several orders of magnitude in com- 
parisons with the corresponding modes for Schwarzschild 
black hole as investigated by BHT. These four modes for 
the electron-positron plasma coalesce with two modes of 
BHT and a single mode of SK and so yielding BHT and 
SK results. Here forlm(k) > corresponds to damp- 
ing and forlm(k) < to growth. This is because the 
convention we have used is e ikz = e i l Rc ( k )+ iIm ( k )] z , The 
limiting horizon values are taken to be n +s — 10 18 cm -3 , 
T +s = 10 10 K, B + = 3x 10 6 G, 7 9 = 4/3, and r) = 0.9 with 
black hole mass M = 5M Q . 



► 




Figure 3: Top: Real part of Alfven damped mode for the 
electron-ion plasma. Bottom: Imaginary part of damped 
mode. 

2. Electron-ion Plasma 

In this case there exist four modes, two of which are 
shown in Fig. 2 and Fig. 3, are growth and damped re- 
spectively and are different modes. The other two modes 
shown in fig. 4 are equivalent to that of Buzzi et al. [6 
and is complex conjugate growth and damping modes. 
The differences in the magnitudes of the oj c x and w C 2 
apparently lead to take the frequencies from their neg- 
ative (and therefore unphysical) values for the electron- 
positron case to positive physical values for the electron- 
ion case. These changes are thus because of the differ- 
ence in mass and density factors as between the positrons 
and ions. These four modes for electron-ion plasma are 
analogous to the four modes obtained by BHT in the 
Schwarzschild case. It is evident that the growth and 
damping rates are independent of the frequency, but de- 
pended only on the distance from the black hole horizon 
through a . The limiting horizon values for electron-ion 
plasma are chosen as n+x = 10 18 cm~ 3 , T + i = 10 10 K, 
n +2 = 10 15 cm- 3 , T+2 = 10 12 K, n = 0.9 , and the mass of 
the black hole M = 5M n . 



Figure 4: Top: Real part of Alfven growth mode for the 
electron-ion plasma. Bottom: Imaginary part of growth 
mode. 

9 High Frequency 
Transverse Modes 

1. Electron-positron Plasma 

There exist four high frequency electromagnetic modes 
for the electron-positron plasma. Fig. 5 shows two modes 
having same real parts equivalent to the real part of Fig.l. 
The mode shown in the center is completely a growth 
mode but the bottom mode is damped for higher fre- 
quency and a Q > 0.5, but growth for a Q < 0.5 and lower 
frequencies. Thus at a distance from the horizon corre- 
sponding to a < 0.5, it appears that energy is no longer 
fed into wave mode by the gravitational field but begins 
to be drained from the waves. The other two modes 
shown in Fig. 6 are not complex conjugate although the 
real part is equal. The center mode shown in Fig. 6 is 
damped for the entire frequency domain but the bottom 
mode is damped for most of the frequency domain and 
for a D > 0.3 and shows growth for lower frequencies and 
for a < 0.3. Comparing the center and bottom modes of 
Figs. 5 and 6 it is quite clear that in the special relativistic 
limit the mode shown in Figs. 5 and 6 are growing and 
damped respectively and coalesce into only one purely 
high frequency mode in the special relativistic case [3 1 j 
for the ultrarelativistic electron-positron plasma. These 
four modes reduce to the three modes of Buzzi et al. [6] 
in the Schwarzschild case. 
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Figure 5: Top: Real part of the two high frequency growth 
and both growth and damping modes for the electron- 
positron plasma. Center: Imaginary part of the growth 
mode. Bottom: Imaginary part of the mode showing both 
damping and growth. 

2. Electron-ion Plasma 

Like the electron-positron plasma, the electron-ion plasma 
admits four high frequency modes. Fig. 7 shows two 
modes as electron-positron plasma. The mode shown 
in the center is a growth mode but the bottom mode 
is growth all the frequency and for a > 0.1, but shows 
growth very closed to the black hole horizon. Thus closed 
to the horizon corresponding to a a < 0.1, it appears that 
energy is no longer fed into wave mode by the gravita- 
tional field but begins to be drained from the waves. The 
third mode shown in the center of Fig. 8 is damped for 
all the frequency domain but the fourth mode shown in 
the bottom of Fig. 8 is damped for most of the frequency 
domain and for a > 0.1 and shows growth for lower fre- 
quencies and for a < 0.1. From these four modes it is 
clear that in the special relativistic limit they coalesce into 



only one purely high frequency mode in the special rela- 
tivistic case investigated by SK[3T] for the ultrarelativis- 
tic electron-positron plasma. The growth and decay rates 
obviously depend on frequency, unlike the corresponding 
Alfven modes. 



[31] where only one purely real mode was found to ex- 
ist for both the Alfven and high frequency electromag- 
netic waves for ultra-relativistic electron-positron plasma. 
In general relativistic limit some new modes arise for 







Figure 6: Top: Real part of the two high frequency damp- 
ing and both damping and growth modes for the electron- 
positron plasma. Center: Imaginary part of the damped 
mode. Bottom: Imaginary part of the damping and 
growth mode. 



10 Concluding Remarks 

In this paper the linearized two-fluid equations are used to 
obtain the dispersion relation for transverse waves closed 
to the event horizon of Schwarzschild Anti-de Sitter black 
hole. In the limit of zero gravity our results reduce to 
those in special relativity obtained by Sakai and Kawata 



Figure 7: Top: Real part of the two high frequency growth 
and both growth and damping modes for the electron-ion 
plasma. Center: Imaginary part of the growth mode. Bot- 
tom: Imaginary part of both growing and damped mode. 

the Alfven and high frequency electromagnetic waves due 
to the black hole's gravitational field. For the electron- 
positron plasma, the damping and growth rates are smaller 
in comparison with the modes for Schwarzschild black 
hole by several orders of magnitude, compared with the 
real components of the wave number. From the vari- 
ous modes for each of the wave types, the damping and 
growth rates are smaller. A change in the radius of the 
black hole, however, does not make any significant differ- 
ence in the result for the field perturbations as it is clear 
from the high frequency waves shown for each fluid com- 



poncnt. For the Alfven waves, the damping and growth 




Figure 8: Top: Real part of the two high frequency damp- 
ing and both damping and growth modes for the electron- 
ion plasma. Center: Imaginary part of the damped mode. 
Bottom: Imaginary part of the damping and growth mode. 

rates are obviously frequency independent, but are depen- 
dent on the radial distance from the horizon as denoted by 
the mean value of the lapse function ao . This is of course 
not the case for the high frequency waves. In that case 
the rate of damping or growth is dependent on both fre- 
quency and radial distance from the horizon. So we make 
conclude that the damping modes demonstrates the en- 
ergy drained from the waves by the gravitational field and 
the growth rates indicate that the gravitational field is, in 
fact, feeding energy into the waves. In the limit £ 2 — » oo 
our study provides the results for the Schwarzschild black 
hole case [6]. 

For a multi-component charged plasma in the pres- 
ence of electromagnetic fields, Marklund et al [32] have 



analyzed an exact 1+3 covariant dynamical fluid equa- 
tions in a Friedmann- Robertson- Walker universe and then 
have considered the matter field of an ion-electron plasma 
with zero average pressure to an Einstein-de Sitter uni- 
verse. The cyclotron damping of the gravitational wave 
produced from resonant interaction with (collisionless) 
plasma particles have also studied by Servin et al |39j and 
Kleidis et al [34 1 135]. while the nonlinear interactions be- 
tween gravitational radiation and modified Alfven modes 
have investigated by Forsberg et al [33] in astrophysical 
dusty plasmas. 

In the context of cosmology, there are also works on 
fluid dynamics and kinetic gas theory such as gas kinetics 
in the Friedmann-Lemaitre- Robertson- Walker (FLRW) mode 
citethirty one, the behavior of matter in the presence of 
electromagnetic fields with plasma effects [3T, 38, i40l |4"T1 
[42], [43] [44] . Therefore, the general relativistic treatment 
of plasmas, both in astrophysics as well as in cosmology, 
seems to be a field open to investigation. 
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